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Abstract 

We propose a characteristic representation of one-dimensional and 
2-state, 3-neighbor cellular automaton rules, which describes an effec- 
tive form of each rule after many time steps. Simulated results of the 
representation show that complex structures of Class IV rules come 
from their aspects of Class II and III. Class IV-like patterns can be 
generated by characteristic functions of a linear combenation of Class 
II and III rule functions. 
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1 Introduction 



Cellular automata (CAs) are discrete dynamical systems with a lattice of 
sites. Each site has a finite set of possible values which evolve synchronously 
in discrete time steps according to identical rules. Wolfram |IJ|| undertook 
valuable works to search systematically for interesting phenomena and pos- 
sible behavior of one-dimensional automata. He found essential four types 
of behavior: homogeneous (Class I), periodic (Class II), chaotic (Class III) 
and complex (Class IV). This classification scheme has been discussed using 
several characteristic parameters [3- 11]. In the previous work[ll] of one of 
the authors and collaborators, characteristic parameters are calculated in all 
rules of the simplest one-dimensional CAs and applied to give quantitative 
grounds for the Wolfram's phenomenological classification. The simplest CAs 
consist of a line of sites with each site carrying a value or 1. The values of a 
particular site is determined by the previous values of its nearest neighbors. 
They are called one- dimensional and 2-state, 3-neighbor CAs or elementary 
CAs. Rules 54 and 110, following Wolfram's numbering^]], are thought to be 
in Class IV||. Typical patterns generated by them show long transients and 
soliton-like phenomena(Fig. 1) . According to the previous paper|Tl|, -Dj//, 



the spreading rate of difference patterns, is the most powerful parameter to 
classify the rules. The calculated results are nearly sufficient for categoriz- 
ing the rules to Class I, II and III. The Class IV rules, however, could not 
be distinguished from Class III rules. Their results of Diff are tentatively 
sorted in Table 1. It is remarkable that rules 54 and 110 exist between Class 
II and III groups. This result supports that the formalism is available for 
characterizing Class IV rules and we can conjecture that their features come 
from interference between their aspects of Class II and III. Simultaneously, 
it represents the confines of discussions on a single characteristic parameter. 
Because rule 41, for example, has the close value of -D;// with those of rules 
54 and 110, we cannot distinguish them. So we present a characteristic repre- 
sentation of each rule, which is given by a linear combination of independent 
rule functions. We can use this formalism to classify CA rules because it 
describes an effective form of each rule function after many time steps. 
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2 Characteristic Representation 



We take x k to denote a value of the /cth site in a one-dimensional Cellular 
Automaton. Each site value is specified as an integer or 1 (2-state). If 
there are N total sites, a configuration can be denoted by 

X = (xx,x 2 , ...,x N ) . (2.1) 

The Hamming distance of X from = (0, . . . , 0) is identical to its weight, 

N 

w(X) = J2 x k- (2.2) 

fc=i 

We define a configuration space {X} N which is a set of all distinct config- 
urations represented by the above form. Time evolution of a site value is 
determined by iteration of the mapping 

T (*+i) f(Jt) Jt) Jt) \ 

= fr(4 ] ), (2-3) 

where / is called a rule function which specifies the CA. The subscript r 
determines a neighborhood of the rule function ((2r + l)-neighbor). In the 
following, we abbreviate the subscript because our discussions are restricted 
to 3-neighbor CAs. A rule function is determined by the succeeding 8 values 

/(0,0 J 0) J /(0,0,1),...,/(1,1,1) ) (2.4) 

or /(0), /(l), • • • , /(7) for short. There are 2 s = 256 rules of the elementary 
CAs and rule numbers are defined by 

E/W'2 1 - (2-5) 

i=0 

Following the notation of the previous paper |IT|, if denotes a rule func- 
tion of the first time step, the function of the £th time step can be obtained 
by 

f {t) (x k ) = fW(fU(...fU(x k ) •••)), (2.6) 
which is a rule for (2t + l)-neighbor of x k : 

f( t) (x k ) = f it) (x k ^...,x k ,...,x k+t ) . (2.7) 
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By the tth rule function we can define a mapping on the configuration space 
{X} N , /<*> : {X} N - {X} N such as 

= /waxi,^,...,^)) 

= (2.8) 

with a proper boundary treatment. In this article, we adopt the periodic 
boundary condition. The image of f®, f^\{X} N ), is a subset of {X} N . 
Since the following is a monotone decreasing sequence of f^\{X} N ys, 

{X} N D fW({X} N ) D fW({X} N ) D • • • D /^({X}*) D • • • , (2-9) 

there exists a limit set lim f^({X\ N ), which will be thought to contain 
characteristics of each CA because its asymptotic behavior is fixed to an 
element or a subset of this set. Classification of CA rules will be carried out 
through arrangement of rule functions by similarity of their limit sets. If 
limit sets of two different rules are equivalent or similar for sufficiently large 
total sites, we can recognize these rule functions to be effectively the same, 
which suggests the existence of a representation of rule functions taking the 
similarity of limit sets into account. 

We take f R ^ to denote the rule function of rule R of the elementary CAs. 
After t time evolutions, the domain of the mapping f { R ] to the (t + l)th 
time step is f^({X} N ). Because : {X} N — > {X} N is a many-to-one 
mapping in general, we define multiplicity of the mapping for an element X® 
of jr({X} n ) as m,R (x^^j. Moreover ur (i,X^^j denotes how many times 
each configuration of three sites, (0, 0, 0), (0, 0, 1), • • • , (1, 1, 1) (i — 0, 1, . . . , 7, 
respectively), appears in a configuration X®. Then the following equations 
are trivial: 

E rn R (X<«>) = 2", E n R (i, X®) = N. (2.10) 

xl*)ef%\x N ) i=0 

Using these parameters, the appearance probability of each configuration at 
the tth time step is obtained by 

P R t) W = ¥ ^ E m fl (*<*>) n fl (<,*<*>). (2.H) 
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Then the average weight of a configuration X^ +1 ^ is given by 



w R 



N 



i=0 



(2.12) 



Modifying the rule function /j^, we define the tth characteristic function 
f R \ The following conditions are imposed on it: 

(i) The mapping f R constructed from f R is defined on {X} N . 

(ii) The value of a characteristic function is limited between and 1: 



o < if (<) < i 



for i = 0,l,...,7. 



(2.13) 



(iii) The average weight of a configuration contained in the image of the 

mapping f^ R is equal to the above average weight of X^ +1 ^. 
The first condition means that the appearance probability Pr (i) is equal to 
|. According to the second, a characteristic function cannot remain in CAs 



anymore. The image of f R is a subset of the unit iV-cube. As we discuss 
in later, a probability interpretation must be given. Although the images 

fR +1 \{X} N ) and f^{{X} N ) are generally not identical, we must define 

f R as the third condition is satisfied. From these conditions, we have 



w R 



N 



o^2-fR 

8 i=0 



(*) 



(2.14) 



Now we propose the following definition of a characteristic function: 

f R (i) for P R (i)>h 



SP R (i)f R (i) + 



j for P R U)>b 



(phCj')-s)/hCj) 



2 for P R (j)>^ 

for P fl (i)<i 



8PrW) 



(2.15) 



When t = 0, is identical with the rule function f R ^ because there is no ef- 
fect of the time evolution. Let us define a function space j /|/ : {^} 3 — > [0,1]}. 
Then the characteristic function and all rules of the elementary CAs, from 
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rule to rule 255, are contained in it. If a product of elements are defined 
by 

(/ • 9) = f[ (/ (0 '9 (<) + (i - / (0) • (i - 9 (0)) - (2-16) 

i=0 

the characteristic function can be expanded by the elementary CA rule func- 
tions as follows: 

255 

/T = E(/r«/S?)/r- (2-17) 

r=0 

The coefficient • means a rate of change from the original rule R to 
rule r after t time evolutions. We call this expanded form as a characteristic 
representation of the rule function. In order to show adequacy of this rep- 
resentation, we calculated the coefficients for all independent even number 
rules of the elementary CAs by computer simulation. Some typical examples 
are listed in Table 2. All Class I rules change to rule 0. Most of Class II 
rules change to itself or typical Class II rules, say 2, 4, 34. Class III rules 
almost stay in the same Class. It is remarkable, however, that the Class IV 
rules 54 and 1 10 have larger changing rate to Class II rules than Class III 
rules (Table 3). These results indicate that their complex structures come 
from the interaction between aspects of Class II and III rules. Repetitions of 
limited patterns as the aspect of Class II are affected by chaotic variations 
as the aspect of Class III, and then complex patterns, soliton-like behaviors 
for example, are generated. In fact, as we will show in the next section, such 
Class IV- like patterns can be generated by characteristic functions of a linear 
combination of Class II and III rule functions. Now we conclude that the 
Class IV rules exist between Class II and III rules, that is, the so-called edge 
of chaos. 

3 Discussion 

Because the tth rule function of the elementary CAs is determined by 
2 2m values, it is difficult to obtain its exact algebraic form after many time 
steps in general. Instead, we considered the characteristic representation 
which took into account the effects of the time evolution not directly but 
through changes of the appearance probabilities of the 3 sites configurations. 
It has an expansion form by all rule functions so that each coefficient means 
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a changing rate from the original rule to the corresponding one. This repre- 
sentation makes clear features of Class IV rules of the elementary CAs as well 
as can be used to classify the CA rules. As noted in Section 2, the calculated 
results showed that the changing rates from Class IV rules to Class II rules 
are larger than the rates from Class III rules to Class II rules. Using the 
characteristic representations, we can generate interference patterns such as 
Class IV rules. If we regard our representation as an extension of CAs to a 
stochastic model, each value of the characteristic function f(x, y, z) can be 
interpreted as the probability that the site y will have a value 1 at the next 
time step. For example, we take such functions that 

f = -^fi + -f90, ^=2/232 + 2/90, (3.18) 

where rules 4, 232 and 90 are typical rules of Class II and III, respectively. 
These functions generate complex patterns as Fig. 2 which are analogous 
to the patterns of Class IV rules. Such interference patterns support the 
conjecture presented in Section 1, that is to say, features of Class IV rules 
comes from interference between the aspects of Class II and III rules. 

Away from the theme of the classification of CAs, this dynamical system 
is very interesting. A pattern construction process can be shown by a con- 
tinuous change of proper elements of f(x, y, z). For example, rule 46 changes 
to rule 110 by varying the value of /(l, 1,0) from to l(Fig. 3). Moreover, 
this formalism will be able to be applied to more complex systems such as 
5- neighbor CAs (now in preparation). 
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